DTP/94-53 
hep-th/M121"91 

December 1994 



Boundary i^-matrix for the quantum 
Mikhailov-Shabat model 

J. D. Kinflfl 

Department of Mathematical Sciences, 
University of Durham, Durham DH1 3LE, U.K. 



ABSTRACT 

We present complete solutions of i^-matrix for the quantum Mikhailov-Shabat 
model. It has been known that there are three diagonal solutions with no free pa- 
rameters, one being trivial identity solution, the others non-trivial. The most general 
solutions which we found consist of three families corresponding to each diagonal so- 
lutions. One family of solutions depends on two arbitrary parameters. If one of the 
parameters vanishes, the other must also vanish so that the solutions reduces to triv- 
ial identity solution. The other two families for each non-trivial diagonal solutions 
have only one arbitrary parameter. 
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I. Introduction 



The Yang-Baxter equation whose solution is called i?-matrix appears frequently in 
physics, eg., factorizable scattering theory on a line and solvable lattice models. 
The actual adaptation of i?-matrix to each of these theories requires some more 
ingredients. 

In the factorizable scattering theory, the solutions of the Yang-Baxter equation 
can be used to compute the iS-matrix, which is given by Z?-matrix multiplied by 
overall factor chosen to satisfy crossing relation and unitarity, etc. Here the spectral 
parameter is related with coupling constant as well as the rapidity parameter of the 
particles and the deformation parameter is related with coupling constant. In the 
solvable lattice theory, the L-operator is given by f?-matrix and the bulk Hamiltonian 
can be computed from transfer matrix which is given by the trace of products of L- 
operators. 

Generalization of factorizable scattering theory to a half-line lead to the discov- 
ery of the reflection equation whose solution is called i^-matrix[|l|]. This equation 
describes interactions related with boundary. By similar way as in the bulk theory, 
boundary ,S-matrix can be computed from K- matrix [fj, |3|. This idea was also used 
to prove the quantum integrability of open spin chains in the framework of quan- 
tum inverse scattering theory[Q. This gives us the Hamiltonian including boundary 
terms. 

At the early stage of development, the solutions of the Yang-Baxter equation were 
computed in a rather direct way, ie., by solving the functional equations component- 
wise. The substantial number of solutions were compiled in Ref. before the general 
algebraic solutions were incidentally obtained^ [7], 

The most general non-diagonal solutions of reflection equation for XXX, XXZ 



and XYZ type i?-matrix which are 4x4 matrices were solved in Ref. [|T], f|, || |10 
The diagonal solutions for bigger size R- matrices were found in Ref. [jTlfl for the 
Zamolodchikov-Fateev model[|T2J, in Ref.|]T3[ for the Mikhailov-Shabat model and in 
Ref. || for A^> model. Spectral parameter independent solutions for various models 
are considered in Ref. |l4j . 
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We consider the most general non-diagonal solutions of reflection equation for 
the 9x9 i?-matrix, especially corresponding to the quantum Mikhailov-Shabat 



model [15]. This model is interesting because it is one of the relativistic quantum field 



theory with single scalar field having non-trivial families of boundary interaction! 16 



This model was known to have two non-trivial diagonal solutions of X-matrix 
as well as trivial identity solution [I3|. These solutions have no free parameters. The 
most general solutions which we found consist of three families corresponding to each 
diagonal solutions. One family of solutions depends on two arbitrary parameters. 
If one of the parameters vanishes, the other must also vanish so that the solutions 
reduces to trivial identity solution. The other two families for each non-trivial di- 
agonal solutions have only one arbitrary parameter and have upper-lower triangular 
structures. 

The plan of this paper is as follows. This introduction is the section I. In section 
II, we consider the whole 9x9 equations resulting from matrix reflection equation 
and describe their general features. In section III, we solve the equations to obtain 
solutions. Finally, we make some discussions in section IV. 
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II. Reflection Equation 

The bulk S'-matrix of integrable quantum field theories can be described by the R- 
matrix which is the solution of the Yang-Baxter equation. The i?-matrix acts on 
V <S> V, where V is (D 3 for the Mikhailov-Shabat model, since this model has 3x3 
Lax-pair representation even though it involves only single scalar field. 

R 12 {u)R 13 {u + v)R 23 {v) = R 23 {v)R 13 {u + v)R 12 (u), (1) 

where Ri 2 (u), R\ 3 (u) and R 23 {u) act on V <8> V ® V, with R\ 2 (u) = R(u) (g> 1, 
R 23 (u) = 1 ® R(u) etc. 

The i?-matrix for the Mikhailov-Shabat or A 2 model was obtained in Ref.fl5 
and it was used to compute S'-matrix in Ref . [|T^] . 
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/(«) 



e u -ll _ g M +»? _|_ e «-3»? _ g -3r; _ g n-5r; _|_ g 5J7 

-e _u + 1 + e" n+4r? - e 4 ", 



-1 + e u + e" 4,? - e u " 4,? . 



This f?-matrix is a meromorphic function of e u . It has only combined PT symmetry. 

V 12 R l2 (u)P l2 = R 12 (u)^\ (4) 

where V%j is the permutation operator on Vi ® defined by "P(x ® y) = (y <E> x) and 
denotes transposition in the i-th. space. 
The boundary versions of the Yang-Baxter equation so-called reflection equations 
for PT symmetric f?-matrix are given by |fl~3| 

R 12 (u - v) K_ (u)R$ 2 (u + v) K_ (v) (5) 
=K- (v)R 12 (u + v) K_ (u)R^ 2 t2 {u - v), 



Ru{-u + v)(K + ) h (u) M~ l R% t2 {-u — v — 2p) M (K + ) t2 (v) 

= (K + ) t2 (v) M Rn(-u -v-2p) AT 1 (K + ) h (u)R\f(-u + v) 



(6) 



1 2 

where K_ (u) = K-(u) <8> 1, K- (u) — 1 ® K-(u), etc. and M is given by V matrix 
defined by the crossing unitarity relation. 



R 12 (u) =V R%{-u -p)V-\ M = V l V = MK 



(2) 

For A 2 } model, p = —677 — m and the matrix V is given by 



V 



In practice, if iT_ (u) is a solution of (|5|) then 

K + (u) = Ki(-u-p)M 

is a solution of (0). 



(7) 



(8) 



(9) 
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Now we write down the 9x9 component equations in the following parametriza- 
tion of X-matrix. 
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1111 
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K(u,rj) 



Unknowns 

zzcd = zzcd 

xxbg + yxed = yybg + xyed 
yybg + xyed = xxbg + yxed 
zzcd = zzcd 
Unknowns 

xzcb + zycg = yzeg + zxeb + xzbd 
zxcd = yzbg + zxbb + xzed 
xzbg + zybb + yzed = zycd 
xzeg + zyeb + yzbd = yzcb + zxcg 
xzcb + zycg = yzeg + zxeb + xzbd 
zxcd = yzbg + zxbb + xzed 
xzbg + zybb + yzed = zycd 
yzcb + zxcg = xzeg + zyeb + yzbd 
Unknowns 

xxce + zzcf = xxce + zzcf 

xybg + xxec + yyef = xybg + xxec + yyef 

yxgb + xygb + xxae + yyae = yxgb + xygb + xxae + yyae 

yxbg + xxe/ + yyec = yxbg + xxe/ + yyec 

yyce + zzc/ = yyce + z^c/ 

Unknowns 

xyea + (3zee + z^/ee + yybg = yxea + z/3ee + xxbg + ^zee 
xyea + (3zee + Z7ee + yybg = yxea + z/3ee + xxbg + 7^ee 
Unknowns 

xzbg + axee + + x(3ec + yzef = xace + (3xcc + zyc/ 

azgg + yx/e + xx</6 + zflfc + + jzff = xyce + (3zcc + zjcf 

xzeg + axbe + (3xeb + x/36c + yzbf = xj3cb + zxcg 1 



(10) 
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1122 : azag + yxge + xxab + z/3gc + (3zgd + jzgf = xxcb + zacg 

1131 : azgg + yxde + xxgb + zfidc + (3zfd + 72c// = z(3cd 

1211 : xz6g + axee + /3x66 + xf3ec + yze/ = xace + (5xcc + zycf 

1223 : xyba + aze& + (3zbe + 276c = zae& + xybd 

1311 : azgg + yxfe + xxg& + zj3fc + (5zdd + jzff = xyce + (3zcc + z7 c / 
1313 : yyfe + xygb + zzfc = yyfe + xygb + zz/c 

1333 : azgg + xyde + ftzdf + yygfr + 7z/d + z^dc = z^cd 

2111 : xzeg + ax&e + (3xeb + x/36c + yz&y = x(5cb + zxcg 

2121 : xxea + zzee + yxbg = xxea + zzee + yxbg 

2132 : yx&a + azeb + z/36e + jzbe = zaeb + yxM 

2211 : azag 1 + yxge + xxab + z(3gc + (3zgd + 7%/ = xxcb + zacg 

2233 : azag + xyge + yyab + /3zg/ + 7zga! + Z7gc = yycb + zacg 

2312 : xyfoa + azeb + /3z6e + z7&e = zaeb + xyM 

2323 : yyea + zzee + xybg = yyea + zzee + xybg 

2333 : yzeg + aybe + xzbf + jyeb + 7/760 = 7/706 + zycg 

3111 : azgg + yxde + xxgft + z(5dc + (3zfd + jzdf = zflcd 

3131 : xxfe + yxgb + zzfc = xxfe + yxgb + zz/c 

3133 : azgg + xg/e + (5zf f + ggg6 + ^zdd + z^fc = yxce + z(5cf + ^zcc 

3221 : yxba + azeb + z(3be + 7z6e = zaeb + yxbd 

3233 : yzbg + ayee + xzef + 7y66 + y7ec = yace + zxcf + 7ycc 

3313 : azgg + xyde + yygb + fizdf + jzfd + zjdc = zjcd 

3322 : azag + xyge + yyab + flzgf + 7zg<i + Z7gc = ggc6 + zacg 

3323 : yzeg + ag&e + xzbf + 7ge6 + g7&c = ^706 + zycg 

3331 : azgg + xyfe + (5zf f + yygb + ^zdd + zjfc = yxce + z/3c/ + 7ZCC 

3332 : yzbg + ayee + xzef + 7y&6 + y~7ec = yace + zxcf + 7ycc 
7 Unknowns 

1213 : aygg + yafe + xagb + zx/c + fiydd + jyff 

= xjbg + ayee + fiybb + xzec + y7e/ 

1231 : aygg + yacfe + xag& + zxdc + /3y/(i + 7yd/ 

= x[5bg + zxbb + y(3ed 

1312 : aygg + yafe + xagb + zxfc + (5ydd + jyff 

= xjbg + ayee + (3ybb + xzec + g7e/ 



7 



1332 : axgg + xade + (3xdf + yagb + 7x/d + zydc 

= yjbg + zybb + X7ed 
2212 : ayag + yage + xaab + zxgc + /3ygd + jygf 

= xaba + axeb + /3x6e + zy&e + yaeg 
3112 : aygg + yade + xag& + zxdc + (3yfd + 7yd/ 

= x[3bg + zxbb + y/?ed 
3132 : axgg + xa/e + (3xff + yag& + ^yxdd + zy/c 

= y/3og + axee + x/3e/ + 7x66 + yzec 
3213 : axgg + xade + /3xd/ + yagb + 7x/d + zydc 

= y^fbg + zybb + x'yed 
3231 : axgg + xafe + (3xff + yagb + 7xdd + zy/c 

= y(3bg + axee + x(3ef + 7x66 + yzec 
8 Unknowns 

1221 : xyeg + aabe + /3ae6 + xxbc + yyft/ 

= xx&a + a(3eb + (3(3be + zzbe + yxey 

1222 : ayag + yage + xaa6 + zxgc + /3ygd + jygf 

= xaba + axe& + (5xbe + zy&e + yaeg 
1321 : axya + yfifb + xzdb + x(3ge + (5xdg + yzge + *yxfg 

= xjeg + aybe + (3yeb + xzbc + yjbf 
1323 : ayga + yzfb + X7<i6 + xzge + f3ydg + y7ye + jyfg 

= y^eg + zyeb + xjbd 
2112: xyeg + aabe + fiaeb + xxbc + yybf 

= xxba + a(3eb + (5j3be + zzbe + yxeg 
2113 : axga + y(3fb + xzdb + x/3ye + /3xdy + yzge + 7x/y 

= X7e<7 + aybe + /3ye6 + xzbc + 2/76/ 
2122 : axaa + y/3y& + xzgb + x/3ae + /3xyg + yzae + 7xyy 

= xaea + ax66 + /3xee + zyee + yafeg 
2131 : axga + y/M) + xzfb + x/3ge + /?x/y + yzge + 7xdg 

= x(3eg + zxeb + y(5bd 
2221 : axaa + y/3y& + xzgb + x/3ae + (3xgg + yzae + 7xyy 

= xaea + ax66 + [5xee + zyee + ya6y 
2223 : ayaa + yzgft + X7y6 + xzae + (3ygg + y7ae + 7yyy 

= yaea + aybb + zxee + xa6y + ^yee 
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2231 : aaga + yxdb + xyfb + xxge + (3a fg + yyge + 7«dg 

= a(3ag + /3/?g/ + yxab + 7/3gd + zzgc 

2232 : axag + lap + (3xgf + yaa& + jxgd + zygc 

= ya&a + ayeb + zx&e + xaeg + 7y6e 
2313 : ayga + + X7(i6 + xzge + /3ydg + y^ge + ^yfg 

= y^eg + zyeb + X7M 
2322 : ayaa + yzg& + x^/gb + xzae + /3ygg + y^ae + jygg 

= yaea + aybb + zxee + xabg + 77/ee 

2331 : ayga + yz<i& + x^fb + xzge + /3y/g + y'yge + 7y<ig 

= y(3eg + ax6e + x/36/ + "fxeb + yz&c 

2332 : yxeg + aa&e + xxbf + 7«e6 + yy&c 

= yyba + a7e6 + zzbe + xyeg + 776c 

3121 : axga + y/3(i6 + xzfb + x/3ge + f3xfg + yzge + 7xa!g 

= x(3eg + zxefr + y(3bd 

3122 : aaja + yxdb + xy/fo + /3a/g + yyge + 7«(ig 

= a(3ag + + yxa& + j(3gd + zzgc 

3123 : ayga + yz<i& + X7/6 + xzge + /3y/g + yjge + ^ydg 

= y(3eg + axbe + x(3bf + 7xe& + yzbc 

3222 : axag + xage + (3xgf + yaa& + 7x5c? + zygc 

= yaba + ayeb + zxbe + xaeg + jybe 

3223 : yxeg + aabe + xx6/ + 'yaeb + yy6c 

= yy&a + a^eb + zzbe + xyeg + 776c 

9 Unknowns 
1322 : aaga + yxfb + xydb + xxge + (3adg + yyge + 70/(7 

= ajag + yyge + xyab + zzgc + fi^gd + 77g/ 
1331 : a/3gg + xxde + /?/3a!/ + yxgb + 7/3/d + zzdc 

= «7gg + yyde + xygb + zzdc + /37/d + 77c?/ 
2213 : aaga + yxfb + xyci6 + xxge + (3adg + yyge + "fafg 

= a^ag + yyge + xya6 + zzgc + (3^/gd + 775/ 
3113 : a(3gg + xxde + (3f3df + yxgb + 7/3/d + zzdc 
= a^gg + yyde + xygb + zzdc + P^/fd + 77^/ 
It is understood that the first, second, third and fourth factors in each terms have 
respectively the arguements u,v,u — v and u + v. Eqs.(1133, 3311) are automatically 
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satisfied. We numbered equations by (ijkl), where ij denotes the rows and kl rep- 
resents the columns of the matrix reflection equation in tensor notation. Equations 
are sorted out according to the number of a, (3, 7, x, y, z, x, y, and z variables. 



Diagonal solutions for K_{u) have been obtained in In this case, it is suffi- 
cient to solve Eqs. (1221, 1331, 1132) only. It turns out that there are three solutions, 
being K_(u) = 1, K_(u) = T + (u) and K_(u) = T~(u), with 



( B 



u 



\ 



G±{u) ) 



'ID 



where 



2 + bf(l 



2 + 6f(l 



2 (-1 iie -3 ") 
1 + e- 6? ? 

2 (1 ± ie 3r ?) 
2-gf{l-e- u Y 01 l + e 6? ? 

These solutions have no free parameters. By the automorphism ( 

for K.(u) follow. 



(12) 
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three solutions 
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III. Solutions 



Now we solve the 9x9 equations listed in the previous section. The basic method 
of solving which we will use is the same as that used when solving the Yang-Baxter 
equation in a direct way||. The known solutions of K- matrix up to now are all 
obtained by this method. The way is, first divide equations by a(u), then take 
the derivative with respect to the variable v and finally set v to zero. This gives 
functional equations involving the variable u only among the elements of i^-matrix, 
which we need to solve. 

We suppose that the solutions to be proportional to identity when the spectral 
parameter u is zero. 

X(0) = X(0) = Y(0) = Y(0) = Z(0) = Z(0) = 0, B(Q) = G(0) = 1. (13) 
We introduce the following functions normalized by a(u) 

X(u) = & Y(u) = M Z (u) = ^M, B( „) = m< ( 14) 
X(u) 

and define their first derivatives at u = as follows. 

x x = X'(u)\ u=0 , yi = Y'(u)\ u=0 , z\ = Z'(u)\ u=0 , (15) 
x 1 = X(u)\ u= o, y~i = Y'(u)\ u=0 , zi = Z (u)\ u=0 , 
bi = B'(u)\ u=0 , gx = G'(u) \ u=0 . 

To begin with, let us suppose that all the elements of if-matrix in Eq.flTUD are 
non-zero and see what happens. Following the procedure described above, we get 
relations between Y(u),X(u) from Eq.(1232) and Y(u),X(u) from Eq.(3212) for the 
2 unknowns. 

vf ^ -xibg + yi ed 

Y{U) = Xl ed-y lb9 X{u) > (16) 
e- r ' +u {x 1 - e u x l + e 3 ^i + e" +u yi) 
e^xi + e u xi - e 3r >yi + e^ +u yi W ' 



x(u) 


Y{u) 


y( 


u) 


Z(«) 




u) 


B(«) 




a(u) ' 


a{ 


uy 


a( 


uy 


ck(m) 


x{u) 


, Y{u) 


v 


(«) 


, Z(«) 
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(u) 
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(«)' 


a(w) 
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- -x x bg + y x ed - 

y u = — -5 -r—X(u), 17 

e- r ' +u {x 1 - e u x x + e 3 ^ + e^yi) -, , 

— — o 5 X(u). 

e 2r >Xi + e u x x - e Zr >yi + e^ +u y x y ' 

Next we go to the equations involving three unknowns. These equations can be used 
to express Z(u) and Z{u) in terms of X{u) and X{u) respectively. Eqs.(1123, 1132, 
1233, 2133) give the same Z(u) and Eqs.(2311, 3211, 3312, 3321) give the same Z(u). 

Z(u) = [_ + J_\ 3v+u _ X(u). (19) 

Until now, the equations involved X, Y, Z and their barred partners seperately. 
However, the equations for 4 unknowns begin to connect the unbarred and barred 
variables. We use these equations to determine X(u) from X(u). It turns out that 
five equations give the same simple relations between X(u) and X(u) 

X(u) = —X(u), (20) 

X\ 

if the following relations holds. 

y~\xx = x x y x . (21) 

Eq . (pi) also give some simplifications to already obtained relations. 

Now we have relations among off-diagonal elements of i^-matrix which can be 
completely determined once we know X(u). To determine X(u) we turn to equations 
involving 6 unknowns, skipping equations involving 5 unknowns. There are sufficient 
number of equations for 6 unknowns to determine X(u). However, the answer is 
rather unexpected in the sense that there is no solutions with every elements of 
.fT-matrix non-vanishing. Let us see how this happens. 

First we notice that from the paired equations like (1121, 2111) we get the 
following relations. 

zt = (22) 

Among 28 equations, we use Eqs.(1211, 2111) to express B(u) in terms of X{u). 
Two equations give two different expressions for the same B(u) function. 

Bl(u) = —r i-Xz x eg - x x be - Yz x bj + Zx x cg + 2X((b, c))) , (23) 

X GO 
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B2(u) = - i( ^_ cc) (-Xzxbg - x x ee - Xb x ec - Y z x e] + Z Vl cf + 2X((e, c))) , 

(24) 

where we introduced Wronskian of two functions as follows. 

((/,<?)) /'// ///'• (25) 

Equating these two different expressions, we can get an expression for X(u). Simi- 
larly we use the equations (2333,3233) to get another expression for X{u) through 
G{u) function. 

Gl(u) = -L- {-Y Zl eg - Vl be - Xz x b] + Z Vl cg + 2Y((b, c)j) , (26) 

G2(u) = — - (-Yzxbg - Vl ee - Y 9l ec - X Zl ef + Zx lC f + 2Y((e, c)j) . 

yi\bb — cc) v 7 

(27) 

So we have two different expressions for the same X{u) function. Requiring these 
to be identical, we get important information about the solutions. There are two 
possibilities. 

xi zi = 0, yi = e~ v x 1 , b 1 x 1 =0, g 1 x 1 = 0. (28) 

4 Ae 2v 

X! Zl = 0, Vl = -e "xi, h = 9i = ^^y- (29) 

It should be noted that the first derivatives of each elements of K-matrix were 
assumed as non-zeros. So the conclusion is, no solution of i^-matrix with every 
components non-vanishing. 

Now we consider the cases when some of the elements of i^-matrix vanish. When 
x~\ is non-zero, above two conditions are still valid so z\ must be zero. When z\ = 0, 
the four equations(1131, 1223, 1333, 2132) imply X(u) = Y(u) = 0. So in this case, 
only X(u),Y(u) and Z(u) can be non-zero. The remaining possibility is only the 
case when x\ itself vanishes. Eq.(3212) for 2 unknowns means Y must also vanish. 
Moreover, we can easily see that Eqs. (2212,1312,3213,3231,3112) for 7 unknowns 
imply either Z = or that X(u) = Y(x) = 0. So in this case, there are two 
possibilities, non-zero X(u),Y(u), Z(u) or non-zero Z,Z(u). 
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In summary, there are three possibilities. Case (I): non-zero Z(u), Z(u), Case 
(II): non-zero X(u),Y(u), Z(u), Case (III): non-zero X(u),Y(u), Z(u). For any 
cases, B(u),G(u) are non-zero. Let me first consider the case (I). 

Case(I): It is easy to see that Z(u) should be proportional to Z(u) eg., from 
Eq.(llll). 

Z(u) = -Z(u). (30) 
z i 

It is useful to start from the Eqs. (21 12,2332) to express the B(u) and G(u) in terms 
of Z{u). They give 

2{{e,b))-z 1 Z{u)be-b 1 eb 

B(U) = 2((e,b)) + b l be (31) 

2-(b 1 + z 1 Z(u))(-l + e- u ) 
2 + 6 1 (l-e u ) 

nl \ 2((e, b)) -Z!Z(u)be-gieb 

G{U) = 2((e,b)) +9l be (32) 

2-( 9l + z x Z{u)) (1 - e") 
2 + ^( e -«_i) 

From Eq.(2123) for 5 unknowns, we get the following relation. 

B(u)z\ee — G{u)z\ee = Z(v)(b\ — gi)ee + 2Z(u)((e, e)). (33) 



Inserting the expressions for B{u) and G(u) in Eqs.(^l|,0) into Eq. (|33|) , we can get 
an expression for Z(v). 

2(-l + e")(l + e")^ 

Z W = TFp » ( 34 ) 

where 

ZF = 4e u + 26i (e u - 2e 2u ) + 2 9l (1 - e u ) + b l9l {\ - 2e u + e 2u ) + z x z x (-1 + 2e u -e 2u ). (35) 

Inserting back this into Eqs. (|31~| , |32"D , we obtain B(u) and G(u) free of Z[u). We 



determined the essential form of the solutions. 

Now we have to see what constraints must be made on the parameters. Only 
one equation is sufficient to determine it. Let us use Eq.(2211) for 6 unknowns. 

z x ag + Zb l9 c - 2Z((g, c)) + Bz l9 d + Gz l9 f = 0. (36) 
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From above equation, we can get 



9l = e 2 ^ z 1 = e 2 ^. (37) 



In all, the solutions for the case (I) is the following. 



Z(u) - (-l + e 2u ) Zl 



Z(u) = 
B{u) = 
G{u) = 



e 2 "(-l + e 2u )h' 



(2e u + bi (e 2r i + e u - e 2u - e 2r ' +u ))z 1 ' 
2e u + 6i(-l + e 2r ? + e u - e 2,?+u ) 
2e u + 61 (e 2j ? + e u - e 2u - e 2, ?+ u ) ' 
e u (2 + - e 2r > - e u + e 2 ^ +u )) 



2e u + &i(e 2r ? + e u - e 2u - e 2, ?+ u ) ' 

bi and zi are the two free parameters. If 61 is zero, g±, z±, Z\ also vanish so that 
the solution reduces to identity. On the other hand, if z\ is zero, so is z\. Then this 
becomes diagonal situation. 

Case(II): In this case, it is quite useful to notice that the Eqs. (1221-2112, 2332- 
3223, 2231-3122, 1322-2213,1331-3113) which are relevant for the diagonal solutions 
remain intact. So the diagonal part for this case is the same as the diagonal solutions. 
Now it is sufficient to find out the X(u),Y(u) and Z{u). We use Eqs.(1112, 1121) 
for 6 unknowns to determine X(u). Eq.(1121) leads to 



Eq.(1112) leads to 



X(u) = ^-^(be + B(u)eb). (39) 



X(u) = ——^——(B(n)(cc - bb) - ee). (40) 
biec + 2{{c, e)) 



To determine Y{u), Eqs. (2333 or 3233) can be used. Eq.(2333) leads to 



Eq.(3233) leads to 



Y(u) = 2(i^) ibE+Giu)eb) - (41) 



y(n)--J^(G(u)(cc-bb)-ee). (42) 
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Above two different expressions for X(u),Y(u) give the same result. 
For b\ = g% — 0, X(u), Y(u) reduce to the following. 

= ( - 1 + e .)( 1 + e . )yi 

For non-trivial diagonal solutions, X(u),Y(u) become 

(±i + e ^)(-l + e")(l + e ")x 1 
X(M) = 2e" (±i + e^) ' (45) 

where ± corresponds to two non-trivial diagonal solutions in Eq. (]TT|) . 

= 2 (±i + ^) • (46) 

To determine the possible conditions on x\ and y±, we use Eq.(1232) for the 2 
unknowns. It gives the following constraints. 

2/i = e-"xi, (47) 

for Eqs.(||,||) and 

yi = iie- 2 ^, (48) 

for Eqs.( f4"5| , |4"E| ) irrespective of the chosen bi,gi in Eq.(|T^). So there are four combi- 
nations to consider, which we denote (++), (H — ), ( — h), ( ), where the first sign 

refers to b\,g\ and the second sign refers to y\. 

Let us turn to Z{u). Eqs.(1123, 1132, 1233, 2133) give the same result for each 
different choices of y\ and b\,g\. For the trivial diagonal solutions, Z(u) becomes 

+ Q (1 + e") (e 2 ^ + e 2 -) z x 

m = 2e2u (1+e2t?) > (49) 

for the yi in Eq . (f47"D . For the non-trivial diagonal solutions, Z(u) becomes 
_ (-1 ± ie*> + e 2 ") (g T ie") (-1 + e") (1 + e») ^ 

ZN " 2e n (±i + e 3, +U) ' ( 50 ) 
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for the (++),( ) choices and 

, v _ (=P + (-1 ± ie^ + e 2r Q je" ± i e") (-1 + e") (1 + e") gi 

W ~ 2e« (±i + en) (±i + e^+ u ) [ ' 

for the (H — ), ( — h) choices. 

We need to determine the possible conditions on z\. We use Eqs.(1223 or 2132). 
For Z(u) in Eq.(^9|), there does not exist consistent solution for z x , ie. Z\,x\ must 
be zero. For the (++), ( ) choices, 

_ (=p + e")(l=F^-e 2 ")a;i 2 

For the (H — ), ( — h) choices, it turns out there is no non-trivial solution for z±, ie. 
zi,x± must vanish. This is the end of the Case(II). The result is the following for 
the sign choices (++), ( ). 

(rpj + e*)(-l ± ie* + e 2r? )(-l + e u ){l + e u )x 1 



X(u) 
Y(u) 



2e u (±i + e 3r >+ u ) 
(1 ± ie*)(-l ± ie 7 * + e 2?7 )(-l + e u ){l + e u )x l 
2e 2r i(±i + e 3r ?+ u ) 



e~ 3r, - u {Ti + e*) (-1 ± ie* + e 2 *f{l - e u ){e* T ie u ){l + e u )xj 
Z[U) ~ 4(±i + e")(±i + e*»+«) ' (53) 



B{u) 
G{u) 



e u (±i + e 3 * +u ) 
e u (e 3 * ±ie u ) 



(±i + e 3r >+ u ) ' 

The solutions which we found satisfies all 9 x 9 equations. 

Case(III): This case can be solved in the same way as in the case (II) since the 
equations for this case is exactly the same. We have only to replace every unbarred 
variables with barred ones. 

( T i + e^X-l ± ie* + e 2 *)(-l + e u )(l + e u )xi 



X(u) 
Y(u) 



2e u (±i + e 3r ? +u ) 
(1 ± ie*)(-l ± ie* + e 2r ?)(-l + e u )(l + e u )x x 



2e 2 *(±i + e 



3ri+u) 



- e- 3n - u {Ti + e*){-l±ie* + e 2,1 ) 2 {l-e u ){e T1 Tie u ){l + e u )x 1 2 

[U) ~ A(±t + e*)(±t + e 3 ^) ' (5 } 
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(e^±te u ) 

B(U) = , r . 

e u (±i + e 3, ?+") 
e"(e 3T >±ie u ) 
(M) ~ (±i + e^+-)- 
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IV. Discussions 



The most general solutions of K-matvix to be determined from reflection equation 
for known i?-matrix are highly desired. Obviously it would be best to obtain them 
in algebraic way. But until now, only small number of solutions are known for small 
size i?-matrices in a direct way as in the present paper. 

We solved the reflection equation for the quantum Mikhailov-Shabat model. This 
model is known to have two non-trivial diagonal solutions as well as trivial diago- 
nal solution. The most general solutions which we found consist of three families 
corresponding to each diagonal solutions. One family of solutions depends on two 
arbitrary parameters. If one of the parameters vanishes, the other must also vanish 
so that the solutions reduces to trivial identity solution. The other two families for 
each non-trivial diagonal solutions have only one arbitrary parameter. 

This solution can be used to compute the Hamiltonian including the boundary 
terms for the corresponding open spin chains along the lines in Ref.Q. It can also 
be used to calculate the boundary S-matrix of the model along the lines in Ref. 
using the i^-matrix and the bulk S'-matrix. The bulk S'-matrix for this model was 



obtained in Ref. 17] from the i?-matrix. 

On the other hand, the quantum Mikhailov-Shabat model has Lagrangian de- 
scription in terms of single scalar field. The classical boundary S-matrix can also 



be computed in the way pursued in Ref.[18|]. It would be nice if we can see whether 
they give the consistent result. 
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